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A long-range "dark force" has recently been proposed to mediate the dark matter (DM) annihi- 
lation. If DM particles are copiously produced at the Large Hadron Collider (LHC), the light dark 
force mediator will also be produced through radiation. We demonstrate how and how precise we 
can utilize this fact to measure the coupling constant of the dark force. The light mediator's mass is 
measured from the "lepton jet" it decays to. In addition, the mass of the DM particle is determined 
using the mr2 technique. Knowing these quantities is critical for calculating the DM relic density. 

PACS numbers: 12.60.Jv, 95.35.+d 



Introduction. Recently, there have been strong indi- 
cations of indirect detections of dark matter (DM) par- 
ticles. In particular, the payload for antimatter matter 
exploration and light-nuclei astrophysics (PAMELA) ex- 
periment has observed a sharp turn-over of positron frac- 
tion in the cosmic rays in the 10-100 GeV range [H, which 
can be naturally explained by DM-DM annihilation to 
electron-positron pairs with a large "boost factor" . In an 
interesting DM scenario the boost factor is obtained 
through the Sommerfeld enhancement effect Q , the con- 
sequence of a long-range attractive force between DM 
particles. A light mediator adark with an O(GeV) mass 
provides the attractive force. Given its small mass, adark 
must couple to the standard model (SM) very weakly 
(hence the force it mediates is given the name "dark 
force" ) , otherwise its effects should have long been seen. 
Accordingly, its coupling to the DM has to be of order 
unity, to obtain the required annihilation rate. The lack 
of anti-proton excess in the PAMELA experiment Q can 
be explained if adark is so light that its decays to hadrons 
arc kinematically forbidden 0], or if it dominantly cou- 
ples to the SM leptons 0. 

The above DM scenario also has unusual signatures at 
the LHC Q. If adark decays within a collider detector 
and dominantly to leptons, we will be able to see "lep- 
ton jets", where two or more leptons are boosted and 
collinear to one another. Once events with lepton jets as 
well as large missing energy are observed, it is critical to 
test whether they come from the same theory that ex- 
plains the positron/electron anomalies. The purpose of 
this letter is to discuss how to extract the relevant infor- 
mation from the LHC measurements to test this new class 
of DM models. In particular, we demonstrate how to 
measure the DM particle mass, the light mediator mass, 
and the coupling constant g of the DM-adark interaction. 
Once these quantities are determined, we will be able to 
calculate DM-DM annihilation rate. This allows us to 
compare with the results of DM searches as well as to 
verify if we obtain the correct DM relic density assuming 
it is thermally produced. 

To be able to perform these measurements, it is crucial 
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FIG. 1: The schematic Feynman diagram of the dark matter 
radiating a light mediator adark, which decays to two leptons. 

that the DM particle can be produced at the LHC, which 
we assume to be the case. Once it is produced, it has a 
significant probability to radiate an extra particle adark 
(Fig.[T]), because adark is light and the DM-adark coupling 
is large. Therefore, for whatever process to produce the 
DM particle, there is a corresponding process with an 
extra adark produced through final state radiation. The 
production rate ratio of the two processes in general de- 
pends on M adark , the coupling g and the DM mass muM, 
which can be written as 

crjpp -> JT DM DM a^k) 
cr(pp -> AT DM DM) 

Here X represents visible particles; C is a process- 
dependent coefficient. This q 2 dependent result is known 
as the Sudakov double logarithm. If other particles 
charged under the dark force arc produced, they can also 
radiate an adark i n which case mDM in Eq. |T]) should be 
replaced by the relevant mass. As we will show shortly, 
M Qdark and todm can be determined independently, which 
enables us to further extract the coupling. 

To be concrete, in this letter we envisage that the 
LHC signal events are produced in a two-step process. 
The first step is entirely analogous to ordinary models 
with missing particles, such as the Minimal Supersym- 
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metric Standard Model (MSSM) with R-parity or Uni- 
versal Extra Dimension with KK-parity. We will adopt 
the nomenclature of the MSSM although a general setup 
is understood. Large cross sections are achieved when 
colored particles are produced and subsequently decay to 
the lightest supersymmetric particle (LSP) in the MSSM. 
Due to the small coupling between the DM and the SM, 
these decays are not affected. In addition, we assume the 
MSSM LSP is heavier than the DM, which is R-parity 
odd. Therefore, in the second step, the MSSM LSP de- 
cays to the DM plus some extra visible particles. In the 
following we exemplify our method by analyzing a simple 
model in the above scenario. After giving the essential 
ingredients of the model, we discuss in turn the mea- 
surements of the mediator mass, the DM mass and the 
coupling. We conclude the letter with a few discussions. 

The measurments. The visible sector in our model 
is the ordinary MSSM with a Bino-like neutralino LSP, 
Xi- The dark sector contains a f/(l)dark gauge group, 
and two Higgs superfields with opposite charges under 
y(l)dark- The dark sector interacts with the ordinary 
MSSM sector through a small kinetic mixing between 
the £/(l)dark gauge superfield and the J7(l)y gauge su- 
perfield. We assume that the dark sector LSP is the 
dark Higgsino, Xdark, which is our DM candidate. Corre- 
spondingly, the dark gauge boson, which is identified as 
the mediator adark i provides an attractive force between 
two Xdark's. ?7(l)dark is broken by the vacuum expecta- 
tion values of the dark Higgs fields (the lightest physical 
Higgs is denoted as /idark); which provide adark a mass of 
0(1 GeV). Due to the kinetic mixing, adark decays to two 
SM leptons, which are collinear with each other and form 
a "dark gauge boson jet" (or a-jet). For the dark Higgs 
^darkj we assume it has a mass M^ dark > 2M adark . There- 
fore it mainly decays to two adark 's which subsequently 
decay to four leptons forming a "dark Higgs jet" (or h- 
jet). Given the assumption that Xdark is lighter than Xii 
X° mainly decays to Xdark plus /i dark or a dark Xi can also 
undergo a three-body decay to x dar k + ^dark + &dark- The 
ratio of the three-body and the two-body decay widths 
can be measured and used to determine the coupling g. 
For this purpose, it is enough to count the number of 
events containing two /i-jets and the number of events 
containing two h-jets plus one a-jet. 

Given the above setup, we choose the dark matter mass 
to be 600 GeV, consistent with the ATIC results. We also 
fix M Qdatk = 1 GeV and M fedark = 3 GeV. The coupling 
constant g is 0.40 to provide the correct DM relic den- 
sity, which determines the boost factor to be 0(100). In 
the MSSM sector, we choose the masses of the gluino, 
squarks (restricted to the first two generations) and the 
MSSM LSP, xi, to be 1200 GeV, 1000 GeV and 700 GeV, 
respectively. The spectrum is chosen such that the gluino 
directly decays to quark plus squark and the squark only 
directly decays to quark plus x?- We generate the par- 
ton level events in the squark/gluino pair production 



channels with Madgraph/Madevents Q for the LHC at 
14 TeV. The total cross section is 0.84 pb. The 2-body 
and 3-body decays for x? are performed with Calchcp 
[H and all other particles, including adark, ^dark and 
the other super particles, are decayed with BRIDGE 
Here, we assume that adark 100% decays to two muons, 
and we will comment on the case with adark decaying to 
electrons later. The parton level events are further pro- 
cessed with PYTHIA [l(| for showering/hadronization, 



and PGS [Hj for detector simulation. 

The lcpton jets are reconstructed as follows: all muons 
are first sorted according to their pr's. Then we choose 
the muon with the highest pt in the list as a "seed" 
for the lepton jet, and add muons within 0.2 rad of the 
seed muon direction to the lepton jet. Used muons are 
removed from the list. The "seed axis" is fixed in this 
procedure. We repeat this procedure until all muons 
are used. Lepton jets with 2 muons are tagged as a- 
jets and lepton jets with 3 or 4 muons are tagged as 
/i-jets. Events containing untagged muons are discarded. 
An Ht > 500 GeV cut is imposed to reduce the SM 
background, where Ht is defined as the scalar sum of 
all objects' pt (including fa) in the events. We list the 
numbers of background and signal events after cuts in 
Table |TJ The SM backgrounds coming mainly from b b 
and ti final states are negligibly small after imposing a 
mass window cut on a-jets. 
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W/Z's 


Signal 


H T > 500 GeV 


1106 k 


1068 k 


17k 


8303 


No. of muons > 4 


168 


1890 


13 


7650 


No. of lepton jets > 2 


70 


36 


< 1 


6523 


|M oje t - 1 GeV| < 0.1 MeV 


1 


1 


< 1 


5970 



TABLE I: The numbers of background and signal events after 
cuts. A 10 fb _1 luminosity and a 14 TeV center of mass energy 
are assumed for the LHC. The final signal events after cuts 
contain 5939 two-body events and 31 2hla three-body events. 
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The mass measure- 



ments of the particles adark and /idark are straightforward 
since that they are given respectively by the invariant 
masses of the a-jet and the /i-jet. The precision of the 
mass measurements is therefore determined by the pre- 
cision of muon momentum measurements. We simulate 
the toroidal LHC apparatus (ATLAS) detector resolution 
with ATLFAST [lj] , which smears both the magnitudes 
and the directions of the muon momenta, and estimate 
the mass measurement error as ~ M/(100\/7V), where 
M denotes M adark or M^ dark , and N is the number of 
corresponding lepton jets. This precision is much better 
than the other measurements, as discussed below. 

Measuring m^o and m^ dark . To determine the masses 
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of Xi and Xdark i we adopt a method based on the vari- 
able rriT2 [13j on events with 2 h- or a-jets. This is an 
inclusive method since that Xi decays to Xdark in both 
decay chains in every event. Alternatively, one could 
consider more specific processes and methods using more 
kinematic information [l4j]. For 10 fb _1 , we obtain 5941 
events with 2 h/ a-jets after the cuts in Table U We use 
A^Xdark an d t° denote some trial masses for Xi an d 
Xdark and examine the number of events that are kine- 
matically consistent with each (/x^ dark fJ>%°)- As discussed 
in [15j, (/^x dark , /Upjo) is consistent with a given event if 
and only if /i^o > WT^^Xdark)' where tot2 is calculated 
using $t and the momenta of the two lepton jets. Using 
this fact, we easily assign the number of consistent events 
for each point on the (Mx da rk> Mx?) mass plane (Fig.[2|). 
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FIG. 2: Contour plot of the number of events consistent with 
the mass pair (jUxdark' ^x?)' Starting from the top contour, 
the number of events descends by 1000 from the adjacent 
contour. Above the top contour, the masses are consistent 
with all 5941 events. The dashed straight line has fixed mass 
difference, which is tangent to the top contour. 



The top contour in Fig. [5] corresponds to the rriT2 end- 
points, TOT2 X (Mxdark)- We note that the m^" contour 
"curls up" on both ends with respect to a straight line 
with constant mass difference and tangent to the contour. 
This is caused by the presence of upstream transverse mo- 
mentum provided by the squark or gluino decay as well 
as the initial state radiation [l6j |. We count the number 
of consistent events along this line, which is maximized 
around /i£j arjl = f)ij darl (Fig. [3]) . We fit the number of 
events distribution with a 5th order polynomial and take 
the maximum position as our estimation of the DM mass. 
The error of this measurement is estimated by repeat- 
ing the procedure for 10 different datasets and comput- 
ing the deviation from the average value, which gives us 
m^ dark = 613 ± 12 GeV (the input value is 600 GeV) 
and m^o — m^ daik = 101.9 ± 0.8 GeV (the input value is 
100 GeV) 

Measuring the coupling constant. As mentioned, wc 
arc interested in the ratio of effective cross sections of 
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FIG. 3: Number of consistent events along the line with fixed 
mass difference and tangent to the m™| x contour. 



2 /idark and 2 ft-dark 1 Odark events. Taking the combinato- 
rial factor into account, this ratio is equal to twice of 
the ratio R = BR(x? -» Xdark ^dark a da rk)/BR(xi -> 
Xdark /idark)- For a fixed ratio of M, ldark /M adark = 3, 
M a dark < m^o - m &ark and m^o - m % 
obtain an approximate formula 



<C m^o, we 
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I log 2 + 4)x + 4(3 log 2 2 + 4 log 2 + 2)] ,(2) 



with x = 2 log [(r 
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)/(^a dark m^o)]. Note 
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FIG. 4: The ratio R as a function of M adark for different 
values of Afh dark /M adark and fixed values m~o = 700 GeV 
and wj dark = 600 GeV. The ratios are proportional to the 
gauge coupling square, which is chosen to be g = 1 here. The 
black dashed line is calculated using the approximate formula 
in Eq. ©. 

when M adark is significantly below 1 GeV, this ratio R 
becomes very large and one should include higher order 
processes with more adark radiated. We also numerically 
show the ratio R as a function of M adark for a few dif- 
ferent values of Mh dark /M adaik in Fig. 2] As can be seen 
from Fig. 2J the approximate formula in Eq. fits the 
numerical results very well, and will be used for error 
analysis of the gauge coupling measurement. The double 
logarithm dependence on r a or M adark is due to the fact 
that the radiated adark tends to be collimatcd with /idark- 
In our model, the annihilation rate for DM to two 
Qdark's or 1 adark plus 1 /idark is given by the S-wave 
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result (for a heavy dark gaugino) 



av 



0.40 



600 GcVJ 



) 2 x 2.31 x 10- 26 cm 3 /s- (3) 



where the required rate 2.31 ± 0.07 x 10 26 cm 3 /s is in- 
ferred [ll from the WMAP result Vt dm h 2 = 0.113 ± 
0.0034 [18[. An 8% increase in the annihilation cross 
section from the Sommerfeld enhancement is taken into 
account in Eq. ([3]). Hence, we have chosen g = 0.40 and 
TO Xdark = 600 GeV to give the correct DM relic density. 
Then R w 0.052 for M 0darlI = 1 GeV. As stated before, 
this ratio can be measured by counting the number of 
events with 2 /i-jcts and the number of events with 2 
h-jets plus 1 a-jet. The latter number is much smaller, 
which dominates the error for the R measurement. For 
10 fb _1 , we expect 220 2hla events before any cuts. How- 
ever, the a-jets tend to be collinear with the h-jets and/or 
contain soft muons that are not registered by the detec- 
tor. This drastically reduces the efficiency for identifying 
2hla events to about 14% [l9| , which has to be taken into 
consideration when calculating the ratio R. The number 
of events is reduced to 31 after lepton-jet reconstruction, 
which results in an error about 18% for the measurement 
of the ratio R. Alternatively, we can treat a lepton jet 
with 5 or 6 muons as collinear /idark + Gdark from Xi three- 
body decay, which then increases the number of events 
to 70 and reduces the error to 12%. 

Having measured the masses and the ratio i?, we deter- 
mine the dark gauge coupling measurement from Eq. @ 
as g = 0.41 ± 0.03. The error for g is obtained by sum- 
ming in quadrature the errors of all parameters it de- 
pends on. Since M adark , Mft, dark and m^a — madark can 
be measured relatively more precisely, the error is domi- 
nated by the errors on R and m^ datk . Using Eq. (3]), we 
conclude that from the LHC measurements with 10 fb _1 , 
the dark matter relic abundance can be calculated to be 
n c d °" zder h 2 = 0.119 ± 0.033. The result is not as precise 
as the measurement from WMAP, but it is certainly very 
important whether they are consistent with each other. 

Discussions and conclusion. For illustration, we 
have assumed that the light mediator only decays to 
muons. The mediator could decay to electrons as well 
resulting in lepton jets containing electrons. To identify 
such lepton jets, we need to consider electrons isolated 
from hadronic activities. Requiring all leptons to be iso- 
lated from hadronic jets by AR > 0.4 and taking the 
electron acceptance into account reduce the efficiency for 
identifying the signal events to around 34%. The error for 
the coupling measurement is then increased by a factor 
of three. Moreover, if the light mediator only decays to 
electrons, we will have to determine the mediator's mass 
from individual electrons' momenta in an "electron-jet" . 
For this purpose, relatively soft electrons (< 10 GeV) 
are favored because they can be sufficiently separated by 
the magnetic field before they hit the electromagnetic 
calorimeter. Moreover, low momentum electrons arc also 



well measured with the tracker. 

In conclusion, we have described how to test a new 
class of DM models with a long-range dark force at the 
LHC. The light mediator decays in the detector to a lep- 
ton jet, which allows us to measure its mass. The DM 
mass is determined using the ttiti technique. Moreover, 
we note that when the DM particle is produced, it has a 
significant rate to radiate an extra light mediator. There- 
fore, we can extract the coupling constant by measuring 
the radiation rate. This technique can be generalized 
to other models with a dark force. For example, the 
decay chain in the dark sector may be longer than the 
one we have considered, involving more particles charged 
under the dark force. In this case, in order to extract 
the coupling constant, one needs to carefully reconstruct 
the decay chain and include all contributions to the dark 
radiation. This and other variations guarantee further 
studies. 
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